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Abstract 



For a class of dynamical systems, called the axiom-A systems, Sinai, Ru- 
elle and Bowen showed the existence of an invariant measure (SRB measure) 
weakly attracting the temporal average of any initial distribution that is ab- 
solutely continuous with respect to the Lebesgue measure. Recently, the SRB 
measures were found to be related to the nonequilibrium stationary state 
distribution functions for thermostated or open systems. Inspite of the im- 
portance of these SRB measures, it is difficult to handle them analytically 
because they are often singular functions. In this article, for three kinds of 
Baker-type maps, the SRB measures are analytically constructed with the aid 
of a functional equation, which was proposed by de Rham in order to deal 
with a class of singular functions. We first briefly review the properties of 
singular functions including those of de Rham. Then, the Baker-type maps 
are described, one of which is non-conservative but time reversible, the second 
has a Cantor-like invariant set, and the third is a model of a simple chemical 
reaction R <-> I <-> P. For the second example, the cases with and without es- 
cape are considered. For the last example, we consider the reaction processes 
in a closed system and in an open system under a flux boundary condition. 
In all cases, we show that the evolution equation of the distribution functions 
partially integrated over the unstable direction is very similar to de Rham's 
functional equation and, employing this analogy, we explicitly construct the 
SRB measures. 
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LEAD PARAGRAPH 



Characterization of nonequilibrium stationary states in terms of dynamical 
ensembles is one of the main questions in statistical mechanics. Recently, the 
so-called Sinai-Ruelle-Bowen (SRB) measures, which had been studied in dy- 
namical systems theory, were found to be related to nonequilibrium stationary 
ensembles for thermostated or open systems. The SRB measures fully describe 
transport properties of the corresponding nonequilibrium stationary states. Also 
they would provide an important insight about the emergence of irreversibility 
in reversible dynamical systems, since they do not have time-reversal invari- 
ance even when the dynamics is reversible. It is therefore illustrative to know 
exact forms of the SRB measures, but it is difficult to handle them exactly be- 
cause they are often singular functions. In this paper, we study three examples 
of Baker-type maps, which illustrate some aspects of the thermostated and/or 
open systems: One is non-conservative but time reversible, the second has a 
Cantor-like invariant set, and the third is a model of a simple chemical reaction 
such as an isomerization R <-> / <-> P. For those maps, we analytically construct 
SRB measures with the aid of a new method, where the weak convergence of 
measures is converted to the strong convergence of partially integrated distri- 
bution functions (PIDFs) and the evolution equations for the PIDFs are solved 
emplying the analogy between them and de Rham's functional equations. 
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I. STATISTICAL MECHANICS AND SRB MEASURES 



One of the main questions in statistical mechanics is to characterize nonequilibrium 
stationary states in terms of dynamical ensembles (cf. e.g., Refs. |]-f|). Recently, for 
thermostated or open systems, stationary nonequilibrium ensembles were found to be related 
to the so-called Sinai- Ruelle-Bowen (SRB) measures @-|T3|, which had been investigated in 
dynamical systems theory [T^-TU]. In the thermostated systems a fictitious damping 

force mimicking a heat reservoir is introduced to avoid an uncontrolled growth of the kinetic 
energy due to an external driving force. The damping force is chosen so as to make the 
dynamics dissipative while it preserves time-reversibility. As a result of the dissipation, 
there exists an attractor of information dimension less than the dimension of phase space 
and the nonequilibrium stationary state is described by an asymptotic measure, which is 
an SRB measure. The SRB measure fully characterizes the transport properties, such as 
the transport law, transport coefficients and their fluctuations. For example, for the driven 
thermostated Lorentz gas ||, the conductivity tensor was calculated, and Ohm's law and 
Einstein's relation were verified by comparing the averaged current with respect to the SRB 
measure to the external electric field. On the other hand, open chaotic Hamiltonian systems 
with a flux boundary condition fl3HT5| admit a nonequilibrium stationary state obeying 
Fick's law that is described by a kind of SRB measure. This measure again characterizes the 
transport properties. Moreover, an open Hamiltonian system with an absorbing boundary 



condition has a fractal repeller, that controls the chaotic scattering [|12| , P0| -Pq] . The unstable 
manifold of the fractal repeller supports a conditionally invariant measure, which provides 
the long time limits of averaged dynamical functions p3| , |24]| . The interrelation between 
the thermostated systems approach and the open systems approach has been discussed by 
Breymann, Tel and Vollmer ||29|| . In this article, we present an analytical construction of 
SRB measures for three examples of Baker-type maps, which illustrate some aspects of 
the thermostated and/or open systems mentioned above. Now we start with the general 
arguments on the SRB measure. 
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The long-term behavior of a dynamical system is characterized by an invariant measure 
fi on an invariant set A, which describes how frequently various parts of A are visited by a 
given orbit x(t) (with t the time). The invariant measure is said to be ergodic if it cannot be 
decomposed into different invariant measures. Such an ergodic invariant measure \i satisfies 
the ergodic theorem |I|-|3|,|TE|] . In case of a map S, it asserts that, for any continuous function 



<p(x), we have 



1 t / (f(x')fi(dx' 



lirn _^ ¥ ,(5*x) = ^ 7 . (x G A\E with fx(E) = 0) (1) 



A 



A dynamical system typically admits uncountably many distinct ergodic measures and not 
all of them are physically observable. One criterion of choosing a physical measure \i is 
that /i describes the time averages of observables on motions with initial data x randomly 



sampled with respect to the Lebesgue measure /iq jTjyiTj]: 



IT = j A ^)»(dz) i xe ^ E with \ 

•' -Z - V ' I n(dx') W 2 ) > and fi a (E) = 0^ w 



t=0 



Sinai, Ruelle and Bowen showed that a class of dynamical systems, called axiom-A systems, 
uniquely admit such a physical invariant measure (SRB measure) |TB|-[TS|]. Thus an SRB 



measure is one for which the ergodic theorem is true for almost every point, x, with respect 
to the Lebesgue measure /i . 

Axiom-A systems are characterized by the hyperbolic structure, i.e., the existence of 
exponentially stable and exponentially unstable directions which intersect transversally with 
each other. In case of bijective differentiable maps (i.e., diffeomorphisms) , the hyperbolic 
structure is defined as follows: Consider a given orbit S l x (t = 0,±1,±2, •••) and small 
deviations 5x of the initial value x. Note that there are m independent possible directions 
of 5x when the phase space dimension is m. Assume that the deviations along m s directions 
decrease more rapidly than an expontial function e~ A * (A > and t > 0) and that the 
deviations along the other m — m s directions grow more rapidly than an exponential function 
e xt (t > and the same A). An invariant set A is said to be hyperbolic when 1) for any 



x G A, the orbit S*x has the m s stable and m — m s unstable directions as explained above, 2) 
the stable and unstable directions depend continuously on x and 3) the stable and unstable 
directions for the point x are mapped by S l to the corresponding directions for the point 
S l x. 

A point x is said to be nonwandering if the orbit S t x returns indefinitely often to any 
neighborhood of its initial point x. If the set Q of all nonwandering points is hyperbolic 
and the set of periodic points is dense in Q, S is called an axiom-A diffeomorphism. In 
particular, if the whole phase space M is hyperbolic, S is called an Anosov diffeomorphism. 
The Arnold cat map is an example of an Anosov diffeomorphism. 

Invariant measures that are smooth along the unstable directions are called SRB mea- 
sures. Sinai, Ruelle and Bowen showed that, for axiom-A diffeomorphisms, the SRB measure 
is the unique physical measure \x describing the time averages (0) of observables of motion 



with initial data x taken at random with respect to the Lebesgue measure \lq [16-18 . For 



more details on axiom-A systems and SRB measures, see Refs. |16|], [|T7|] and ||19|| . 

In Gibbs' picture of statistical mechanics, a macroscopic state for an isolated system 
is described by a phase-space distribution function, and a macroscopic observable by an 
averaged dynamical function with respect to the distribution. Suppose that the dynamics 
satisfies the mixing condition with respect to the microcanonical distribution /i mc : 



/ (p(x)/j mc (dx) 

lim / ip(S t x)p (x)p mc (dx) = J -^ L 1 , (3) 

c ^ +ooJm / vUdx) 

JM 

where M stands for the whole phase space, <p(x) is a continuous dynamical function and 
Po(x) is a normalized initial distribution function. Then, the system exhibits time evolution 
as expected from statistical thermodynamics. Namely, the distribution function weakly 
approaches an equilibrium microcanonical distribution and the averaged dynamical functions 
approach well defined equilibrium values From this point of view, instead of Eq. (0), it 
is enough to consider the following condition as a criterion of choosing a physical measure 
//: 
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(p(x)fi(dx) 

lim I (f(S t x)p (x)no(dx) = , (4) 

fi(dx) 



A 



where /xq stands for the Lebesgue measure, M is the whole phase space, A the attractor and 
<p(x) and po(x) are respectively a continuous dynamical function and a normalized initial 
distribution function. Sinai, Ruelle and Bowen JIB] showed that the SRB measures for the 



axiom-A systems satisfy Eq. (Q) as well. Hence, the measure satisfying Eq. (Q) will be 
refered to as a physical measure. Since the left hand side of Eq. (f|) is the average (ip)t of tp 
at time t, Eq. (|j) can be generalized to define a physical measure /i for systems with escape 



23,24 



<p(S x) po(x) po(dx) / ip(x)fi(dx) 

lim>)t = t lim. "T* = ^7 > ( 5 ) 

Xm(5'^)po(^)a*o(^) / M<^) 



t— >+oo t^+oo 



M J A' 



where xm stands for the characteristic function of the phase space M and A' is the support 
of \x. The denominator of the middle expression is necessary as the total probability is not 
preserved. When there is no escape, the physical measure defined by Eq. (Q) is supported 
by the attractor and, thus, is invariant. On the other hand, when there is escape, the 
physical measure defined by Eq. (|5|) is, in general, not an invariant measure. However, 
since the ratio f A , <p(x)p,(dx)/ f A ,p(dx) does not evolve in time, such a measure is called 
conditionally invariant f2~3|j2"4|| . We remark that, when they are of axiom-A type, the systems 
with escape also possess "natural" invariant measures supported by the fractal repeller, 
which are specified e.g., by a variational principle []I|, |12| , p2| -p8| . It is the natural invariant 
measure that characterizes the ergodic properties such as the Lyapunov exponents, but not 
the conditionally invariant physical measure defined by Eq. (|5|). 

Now we revisit thermostated and open systems. As explained before, a thermostated 
system is dissipative because of the fictitious damping force. Then, a nonequilibrium sta- 
tionary state is described by an asymptotic SRB measure defined by Eq. (|2|) ]^T2[|. For 
open chaotic Hamiltonian systems with a flux boundary condition, a nonequilibrium sta- 
tionary state obeying Fick's law is described by a measure with a fractal structure along the 



contracting direction [|13|-|T5|1. Since the measure is smooth along the expanding direction 



and can be defined by Eq. (^), it is an SRB measure [|X3j] . In both cases, those SRB measures 
fully characterize the transport properties. However, it should be noticed that the two cases 
are different because the invariant set of an open system is a fractal repeller which is fractal 
in both the stable and unstable directions and, thus, does not support an SRB measure, 
while the invariant set of a thermostated system is an attractor which does support an SRB 
measure. 

By applying Eq. @ for an initial constant density, one obtains a method of constructing 
an SRB measure for a map f30| : 1) Approximate the measure by iterating an initial constant 
density finite times, 2) calculate the average with its result and 3) take the limit of 
infinite iterations. Several methods are also proposed where unstable periodic orbits or 
trajectory segments are used to write down an SRB measure and averages with respect to 



it (cf. Refs. P,|7|,|8] J10| , |2T|p2] , |31|| and references therein). However, it is not easy to evaluate 



the convergence rate of the limits in Eqs. (^))-(0), particularly for non-expanding maps, 
and to extract exponentially decaying terms from an averaged dynamical function {<p)t, 
which are the Pollicott-Ruelle resonances |32|-^,0,^-|5S| . One of the reasons is that the 



long-time limit of the measure can be defined only via the ensemble average as shown 
in Eqs. (§) and (||). For non-expanding maps, the distribution function itself does not 
have a well-defined long-time limit. In other words, an asymptotic SRB measure is the 
weak limit of an initial measure. In an analytical construction of SRB measures which we 
shall explain for three Baker-type maps, the weak convergence of measures is converted to 
the usual convergence (technically speaking, the strong convergence) of partially integrated 
distribution functions. Contrary to the evolution equation of a distribution function (the 
Frobenius- Perron equation), the evolution equation of a partially integrated distribution 
is contractive and possesses a well-defined long-time limit. This equation is similar to de 
Rham's functional equation which was introduced to deal with singular functions such 
as continuous functions with zero derivatives almost everywhere, and its contraction rate 
gives the convergence rate of the averaged dynamical function (<p)t- In Sec. II, we review 
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some properties of singular functions including those of de Rham's functional equation. In 
Sec. Ill, an SRB measure is constructed for a non- conservative reversible Baker map with 
the aid of de Rham's equation. The model illustrates the two fundamental features of the 
thermostated systems, namely the phase space contraction and time reversibility, and we 
discuss the interrelation between the two features. In Sec. IV, a Baker map with a Cantor- 
like invariant set is studied. When there is no escape, the map possesses a strange attractor, 
which is a direct product of the unit interval and a Cantor set. On the other hand, when 
there is escape, the map has an invariant set, which is the direct product of two Cantor sets, 
and is a simple example of an open system with an absorbing boundary condition. Physical 
measures for the map with and without escape are constructed with the aid of de Rham's 
equation and the natural invariant measure for the map with escape is derived. In Sec. V, 
we investigate the properties of a Baker map with a flux boundary condition, which mimics 
a chemical- reaction dynamics with a flux boundary condition (cf. Ref. |[40||). We show that 
an SRB-type stationary distribution describes the reaction dynamics and that the slowest 
relaxation to it is characterized by a decay mode (i.e., the Pollicott-Ruelle resonance), which 
is a conditionally invariant measure. Sec. VI is devoted to concluding remarks. Technical 
details of the construction of measures are presented in Appendixes. 



II. SINGULAR FUNCTIONS AND DE RHAM EQUATION 

Basic tools of our analytical construction of SRB measures are singular functions and 
de Rham's functional equation. Singular functions such as the Weierstrass function or the 
Takagi function were originally introduced as pathological counter examples to the intuitive 
picture of functions. These singular functions play an important role in chaotic dynamics. 
A step towards the analytical treatment of singular functions was given by de Rham p9| . 
who showed many of them can be characterized as a unique fixed point of a contraction 
mapping in a space of functions. In this section, we briefly review the properties of some 
singular functions and the relation between them and chaotic dynamical systems (cf. Refs. 
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47[ and glj). 

The first example of a nowhere differentiable continuous function was given by Weier- 



strass in 1872 |TJ (Fig. |T]): 



W a Jx) 



COs(b n 7TX 



(6) 



n=0 



where 0<a<l,6isa positive number and ab > 1. Moser |42j used the Weierstrass 
function to construct a nowhere differentiable torus for an analytic Anosov system. Yamaguti 



and Hata [f43[] used it as a generating function of orbits for the logistic map and discussed 
some generalizations. Also Weierstrass functions are eigenfunctions of the Frobenius-Perron 



operator for the Renyi map Sx = rx (mod 1) with r a positive integer |35|,44]|. 

In 1903, Takagi gave a simpler example of a nowhere differentiable continuous function 
p[ (Fig. |): 



OO 1 



(7) 



where ip(x) = \x — [x + 1/2] | and [y] stands for the maximum integer which does not exceed 
y. In 1930, van der Waerden gave a similar function which is obtained from Eq. (0) by 
replacing 2" to 10 n ||46||. In 1957, the Takagi function was rediscovered by de Rham 



Some generalizations of the Takagi function were discussed by Hata and Yamaguti |j43lf47 



The function T(x) — 1/2 is the eigenfunction of the Frobenius-Perron operator for the map 
Sx = 2x (mod 1) with eigenvalue 1/2. Also, the Takagi function and related functions were 
found to describe the nonequilibrium stationary state obeying Fick's law for the multi-Baker 
map 
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In the theory of the Lebesgue integral, there appear nonconstant functions with zero 
derivatives almost everywhere, which are sometimes referred to as Lebesgue's singular func- 
tions |48| , |49"|j . One typical example is the Cantor function (Fig. [3]). A more interesting 
example f a (x) (0 < a < 1, a ^ 1/2) is the unique function satisfying 

af a (2x) , x G [0, 1/2] 

(8) 

(1 - a)f a (2x - 1) + a , xe [1/2,1] 
10 



fa(x) 



which is strictly increasing and continuous, but has zero derivatives almost everywhere with 
respect to the Lebesgue measure [^] (Fig. |). Note that such functions do not satisfy the 
fundamental theorem of calculus [[48| : 



fa{l) ~ /«(0) = 1 ^ = ^ f' a {x)dx. 

Jo 

The function f a (x) with real a represents a cumulative distribution function of an ergodic 



measure for the dyadic map Sx = 2x (mod 1) (cf. examples given on p. 626 of Ref. |16 
and on p. 36 of Ref. |pOf| ). The eigenfunctionals of the Frobenius-Perron operator for the 
multi-Bernoulli map and the multi-Baker map can be represented as the Riemann-Stieltjes 



integrals with respect to f a (x) with complex a |36|.|37]]. In Ref. |3q| , it was shown that, for 
a class of piecewise linear maps, the left eigenfunctionals of the Frobenius-Perron operator 
admit a representation in terms of singular functions similar to f a (x). As pointed out by 



Hata and Yamaguti [[43|,|^J], f a (x) is analytic with respect to the parameter a though it 
is a singular function of x, and there exists an interesting relation between the parameter 
derivative of f a (x) and the Takagi function: 



iL Ux) 



= 2T(x) . (9) 

a=l/2 



In 1957, de Rham found that the Weierstrass function, the Takagi function and 
Lebesgue's singular function as well as other singular functions are fully characterized as 
a unique solution / of a functional equation 

f(x) = Ff(x) + g{x) , (10) 

where g(x) is a given function and JF is a contraction mapping with JFO = defined in 
the space of bounded functions on the unit interval [0,1] [|39|. Then, he generalized the 
functional equation (|K]) to describe fractal continuous curves such as the ones by Koch or 
Levy. The contraction mapping means that the inequality 



||^/i-^/2|| <A||/ x -/ 2 
11 



holds for some constant < A < 1 and for any functions f\ and f% } where the function norm 
||-|| is defined as the supremum: ||/|| = sup[ 0il ] The existence of a unique solution of 

de Rham's functional equation (|T0| ) immediately follows from Banach's contraction mapping 



theorem |^9[ and the fact that the mapping / — > T f + g is contraction. Let the mapping 

Fa,p be 

f af(2x) , x G [0, 1/2] 

Faj>f(x) = , (11) 

(Pf(2x-1), xe (1/2,1] 
then the Weierstrass function 1^,2, the Takagi function T(x) and Lebesgue's singular func- 
tion f a (x) are the unique solution of ([TDD respectively for a — ^ — a, g(x) = cos7ra;; for 



a — (5 — 1/2, g(a;) = |sc — [x + 1/2] |; and for j3 — 1 — a, g(x) = a9(x — 1/2) with 9 the step 
function |39[ . For more information on singular functions, see Refs. [[|7]] and |5l] . 



Before closing this section, we remark on Riemann-Stieltjes integrals with respect to 
singular functions, which will appear in the next section. Suppose f(x) is of bounded 
variation and <f(x) is continuous. Then, the Riemann-Stieltjes integral of (p with respect to 
/ is defined by the limit 

r b n 

/ ^x)df{x)= lim E^H/W-ZW}, (12) 



where {x^} is a partition of [a, 6]: a = x < X\ < ■ ■ ■ < x n = b, and Xk-i < £,k < %k 
Since the formula for the integration by parts 



\{x)df{x) + f b f(x)d V (x) = <p{b)f{b)-<p{a)f(a) , 

J a 

holds in general, the Riemann-Stieltjes integral of a function of bounded variation with 
respect to a continuous function can be defined as above. At first sight, the evaluation 



of the Riemann-Stieltjes integral ([12]) seems to be difficult. But, for a class of functions 
obeying de Rham's equation, it is not the case. As an example, consider the Fourier-Laplace 
transformation of f a (x): 



1{j] )= I* e^df a (x) 



12 



By dividing the integral into the ones over [0,1/2] and [1/2,1] and changing the variable x 
to x/2 and (x + l)/2 respectively, we have the recursion relation 



1 rl 

»W 2 Jf Ct/9U / Jv(x+i)/2 f 



I( v ) = / eWdMx/2) + / e^^df a ((x + l)/2) 
Jo JO 

= a [ e ivx/2 df a (x) + (1 - a)e ir > /2 [ e ivx/2 df a (x) 
jo jo 

= {a + (l-a)e^ 2 }I( V /2) , 

where the de Rham equation (|T0|) is used in the second equality. Because 1(0) = Jq df a (x) 
/„(!) = 1, the above recursion relation gives [^ , ^ , pl| , p2]] 



/ e^dUx) = = ]J{a + (1 - a)e^ 2n } . (13) 

Jo n=l 

Note that, for a ^ 1/2, I(2 m n) = I(2n)(^ 0) (m = 0, 1, • • •) and, hence, the Fourier-Laplace 
transformation of f Q does not satisfy the Riemann-Lebesgue lemma: lim^oo I (if) ^ 0, that 
again implies the singular nature of f a . 

Further we notice that the formula ( |T3"D and the Hata-Yamaguti relation (|S[) relate the 
Lebesgue's singular function and the Takagi function to the Weierstrass functions: 

r z \ Im/(27rs) Re/(27rs) -1 ~ , , 
f a (x) = a - 22 W 1/2:2 (2sx) + \J i '- W 1/2 , 2 (2sx) , 

s>0:odd nS s>0:odd 7TS 

T ^) = \~ E 4~ 2 W W (2sx) , 

Z s>0:odd " 6 

where the sums runs over positive odd integers, Wx/ 2>2 (x) is the Weierstrass function (|J) and 
Wi/2, 2 (x) is a singular function obtained from (|6|) by replacing cos to sin. 



III. SRB MEASURE FOR A NON-CONSERVATIVE REVERSIBLE BAKER MAP 

In thermostated systems, dynamics is non- conservative due to the damping force mimick- 
ing a heat reservoir and thus, the forward time evolution is different from the backward time 
evolution. However, it is time reversible |4]-|T2"fl. It is therefore interesting to see how these 
two features are compatible, that we shall study with a simple map. One of the simplest 
non- conservative systems which have time reversal symmtery is given by (Fig. [5]) 
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(x/l,r y) , x G [0,/] 

*{x,y) = { (14) 
((a; - l)/r, I y + r) , x e (1,1] 

where / + r = 1 and < I < 1. The map is non-conservative since its Jacobian takes r/7 

for x G [0,1] and //r for x G (/, 1], both of which are different from I. But the map has a 

time reversal symmetry represented by an involution I(x,y) = (1 — y, 1 — x): = < 1 ) ~ 1 . 

The Frobenius-Perron equation governing the time evolution of distribution functions (with 

respect to the Lebesgue measure) is given by 



p t+ i{x,y) = Up t (x,y) = J dx'dy'5[(x,y) - §(x' ,y')]p t (x' ,y') 

(15) 



f Pt(l x,y/r) , y G [0,r] 



j p t (r x + l,(y-r)/l) , y G (r, 1] 

where U stands for the Frobenius-Perron operator defined by the second equality and S[-] 
is the two-dimensional delta function. Since the map $ is not conservative, the numerical 
factors r/l and l/r different from 1 appear in the last expression. 

A. An SRB measure for the forward time evolution 

First we explain our algorithm to construct SRB measures and apply it to the forward 
time evolution. Our goal is to show that an expectation value of the dynamical function 
ip(x, y) with respect to p t converges for t — > +oo to an expectation with respect to a singular 
measure given below, when the initial distribution function p is continuously differentiable 
in x, and a dynamical function ip is continuously differentiable in y and continuous in x. We 
remark that the convergence rate is controlled by the smoothness of the initial distribution 
function and the dynamical function, and the condition given above is sufficient for the 
exponential convergence. 

The first step in our explicit construction of the singular measure is to express the 
expectation value by the partially integrated distribution function, i.e., 



J[o,i} 2 J[o,i} 2 
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<f(x, 1) P t (x, l)dx - f d y ip(x, y) P t (x, y)dxdy , (16) 

"'[0,l] 2 

where P t (x,y) = f$ dy'p t (x,y') is the partially integrated distribution function, the symbol 
d y stands for the Riemann-Stieltjes integral of P t only with respect to the variable y [|53| 
and the integration by parts is used in the second equality. The evolution equation for P t 
can be obtained easily from Eq. (|T5|): 



I P t {l x,y/r) , ye [0,r] 

Pt + i(x,y) = { (17) 
r P t (r x + l,{y- r)/l) + I P t (l x, 1) . ye (r, 1] 

Partial integration of the distribution function changes the prefactors from r/l and l/r 

respectively to r and I, which are strictly less than unity. Because of this, the evolution 

equation ([17|) is similar to de Rham's functional equation ([T0|) . 

The next step in our construction of the singular measure is to calculate the long time 

limit of P t . Putting y — 1 in Eq. (|17|), we obtain 



P t+1 (x, 1) = r P t {r x + l,l) + l P t (l x, 1) . (18) 

Note that this is nothing but the Frobenius-Perron equation for a one-dimensional chaotic 
map (strictly speaking, an exact map cf. Sx = x/l (for x e [0,/]) and Sx = (x — l)/r 
(for x e (^,1]), which admits the Lebesgue measure as the invariant measure. Therefore, 
the normalization integral Jq 1 dxP t (x, 1) is invariant: 

/ dxP t (x, 1) = / dxP t -i(x, 1) = • • • = / dxP (x, 1) , 
Jo Jo Jo 

and is equal to the long time limit of the partially integrated distribution function Pt{x, 1) 
EjjJ . As will be shown in Appendix A, the convergence rate is A. 

f 1 dx'P (x', 1) + O(A') = / p (x', y')dx'dy' + O(A') . (19) 

JO J\0A} 2 



In order to proceed with the calculation, one needs the following lemma. (For its proof, 
see Appendix A.) 
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Lemma Let T be a linear contraction mapping with the contraction constant < A < 1 . 

And let g t = g^ + ^g^+gf^ be a given function where v is a constant satisfying A < \v\ < 1, 
and g^ = O(A'). Then, the solution of the functional equation 

ft + x = Fft + 9t, (20) 

is given by 

ft = /£> + SfW + O(tA') , (21) 
where fgp and f^' are the unique solutions of the following fixed point equations 

Jt ] = Ff£> + 9 {0) , (22) 

/<? = ^/£> + ^ • (23) 



Now we go back to the equation (|17D for Pt(x, y), which can be rewritten as 

P t+1 = FP t + g t , (24) 

where the contraction mapping T and a function g t are respectively given by 

(lP{lx,y/r), ye[0,r] 
FP(x,y) = \ (25) 
[ r P(r x + l,(y — r)/T) . y G (r, 1] 

and 

!0 y G [0 r) 

= 9 {0) (y) [ p^x'^dx'dy' + OiX 1 ) , (26) 
lP t {lx,l), ye(r,l] VH 2 

with 

f0, ye[0,r] 

g {0 \y) = (27) 

[l . y G (r, 1] 

The contraction constant of the mapping T is A = max(/,r) and Eq. flTP| ) is used to derive 
the left-hand side of Eq. (|2T)|). As the equations (^), fl2"5]), ( P^| ) and fl2"T| ) satisfy the condition 
of the lemma and the contraction mapping T given by (^) is linear, the lemma implies 

16 



P t (x,y) = F t (y) I p (x',y')dx'dy' + 0(tA*) , (28) 
J [0,1] 2 

where Fi(y) is the unique solution of de Rham's functional equation 

.lFt(y/r), ye[0,r] 
F l (y) = fF l (y)+g(°\y) = { (29) 

rFi{{y-r)/l) + l. y G (r, 1] 



By substituting Eq. (p8|) into Eq. ( [lq) and employing the integration by parts, we have 



(p(x,y)p t (x,y)dxdy = (p(x,y) dxdF t (y) / p (x' ,y')dx'dy' + 0(t\ ) . (30) 

[0,1] 2 "'[O,!] 2 "'[O,!] 2 

We remind the reader that Eq. ( [3~0|) holds for any continuous function (p(x,y) and any 
integrable function p (x,y) which are continuously differentiate respectively in y and x. If 
Po is normalized with respect to the Lebesgue measure, Eq. ( |3~0"D gives 



lim / (f(x,y)p t (x,y)dxdy = / ip(x,y) dxdF t (y) . (31) 

t-*ooj[0,l] 2 "'[O.l] 2 

This shows that the physical measure yU p h of the system is given by 

p ph ([0,x) x [0,i/)) =zF,(s/) . (32) 
Clearly it is absolutely continuous with respect to the Lebesgue measure along the ex- 



panding x-direction and, thus, is an SRB measure. As studied in Ref. |pq| , the func- 
tion Fi is non-decreasing, has zero derivatives almost everywhere except for r = I = 1/2 
with respect to the Lebesgue measure and is Holder continuous with exponent 5 = 
— lnmax(/, r)j lnmin(/ _1 , r -1 ) = 1 (i.e., \F\[x) — Fi(y)\ < A\x — y\). The graph of Fi is 
a fractal (Fig. 6), but its fractal dimension is D = 1 as a result of the Besicovich-Ursell 
inequality f56|: 1<D<2 — 5=1. Moreover, the one-dimensional measure defined by Fi is 
a multifractal two-scale Cantor measure, the dimension spectrum D q (— oo < q < +oo) of 



which is given as the solution of |38,55 



There exists an interesting relation between Fi and Lebesgue's singular function f a 
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F l (y) = f l of r -\y) , 

which immediately follows from the fact that the right-hand side obeys the same functional 
equation as i 7 ). Note that, since f r is continuous and strictly increasing, the inverse / r _1 
exisits and is also strictly increasing. As a composite function of two strictly increasing 
functions fi and f' 1 , Fi is also strictly increasing. Because of these singular properties of 
Fi, the physical measure // p i, is singular along the contracting y-direction. 

The physical measure /x p h is mixing with respect to the map $. Indeed, by considering 
P t (x,y) = $ p t (x,y')dFi(y') instead of P t (x,y) = $ p t (x,y')dy' and following exactly the 
same arguments as above, one obtains 

Inn/ *x,v)«fa«F,to), (33) 

*^°°J[0,1] 2 "'[0,1] 2 

provided that (p(x, y) and po(x, y) are continuously differentiable respectively in y and x and 
that po is normalized: / po(x, y)dxdy = 1. By using this fact, the Lyapunov exponents can be 
analytically calculated as follows: The Jacobian matrix for $ is diagonal, and the logarithm 
of the component along the expanding x direction, the local expanding rate A x (x,y), is 



K(x,y) 



- In I , x G [0, 1] 

— lnr . x E (1,1] 

The Lyapunov exponent along the x direction (the positive Lyapunov exponent) is defined 
as the average of A x (x,y) over an orbit starting from some initial point (x,y). Since the 
system is ergodic, the Lyapunov exponent can be obtained from Eq. @ using the measure 
P P h ■ 

j 

A x (^(x,y))= , 

'[o,i 



A 



l T r 

x (<S>,p ph ) = T \im^-J2A x {^ t (x,y))= J^^A x (x,y)dxdF l (y) = -ZlnZ-rlnr . (34) 



Similarly, the Lyapunov exponent along the y direction (the negative Lyapunov exponent) 
is 

\ y (®,p ph ) — liar + r\al . (35) 



18 



The sum of the two Lyapunov exponents is negative 



Aa:($, /iph) + A y ($, /i ph ) = (r - /) In ( - ] < . 



Hence, areas are contracted on average by the map $ and the map possesses an attractor A. 
From Eq. fl3~T|), one finds that the attractor A is the support of the SRB measure /i p h- When 
I 1/2, the two-dimensional Lebesgue measure of the attractor A is zero since the function 
F\ has zero derivatives almost everywhere with respect to the Lebesgue measure. Moreover, 
according to Young's formula |57| (which is the Kaplan- Yorke formula []T6| , pLTL 30 , 5 8 1 for two- 
dimensional ergodic systems), the information dimension of A is given by : 



Az($,/iph) 



= 1 + 



/ In I + r In 



I In r + r In I 



< 2 



|A y ($,/x P h)| 

which is less than two. Therefore, A is a fractal set. On the other hand, since the function 
Fi is strictly increasing, for any non-empty rectangle [x ,x + e) x [yo,Uo + e') (e > and 
e' > 0), we have 



fi p h([xo,x + e) x [y ,yo + e')J = eji^yo + e') - F { (j/ )} > , 

which implies that A fl [x ,x + e) x [y ,Z/o + e') 7^ and, hence, the attractor A is dense 
in the whole phase space [0, l) 2 . This phase-space structure is in contrast to the one of a 
dissipative system usually studied ]30| (see also the next section), where an attractor is a 
Cantor-like set. 



B. An SRB measure for the backward time evolution 

Now we consider the backward time evolution. In the same way as the forward evolution, 
we find that another partially integrated distribution function P t (x,y) = f$ dx' p t (x' , y) 
converges for t — > — 00 and we have 

/ ip(x,y)p t (x,y)dxdy = <f(x,y) dF r (x)dy / p (x', y')dx'dy' + 0(|t| A 1 * 1 ) , (36) 
J[o,i] 2 Ao,i] 2 "'[0,l] 2 
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where t = 0, —1, —2, • • • and a singular function F r is given by F r (x) = 1 — Fi(l — x). As 
before, Eq. fl36|) implies that the asymptotic physical measure /2 p h is given by /2 p h([0,x) x 
[0, i/)) = F r (x)y. The measure /2 p h is then absolutely continuous with respect to the Lebesgue 
measure along y direction and singular along x direction. This correponds to the fact 
that the expanding and contracting directions are interchanged for the backward motion. 
Actually, the measure /2 p h is precisely the one obtained from ^ p h via the time reversal 
operation I: /i p h = //iph- The measure /i p h is again mixing with respect to the backward 
time evolution $~* (t = 0, 1, • • •). And the Lyapunov exponents are calculated as the phase 
space averages of the local scaling rates for the inverse map $ _1 with respect to /2 p h. For 
example, the positive Lyapunov exponent is the /i p h-average of the local expanding rate: 
A y (x, y) = —6(r — y) lnr — 6(y — r) In / and is equal to that for the original map $. 

A y ($ _1 ,/2 ph ) = / Ay(x,y)dF r (x)dy = -rlnr - Ihxl = X x (^,ii ph ) . (37) 

J [0,1] 2 

The negative Lyapunov exponents of the two maps are also the same: 

A x ($-\// ph ) = A^,//^) . (38) 

The equality of Lyapunov exponents for ($,/i p h) and ($ _1 , ft p h) is a general consequence of 
the time reversal symmetry of the system. 

We notice that the natural invariant measure /i p h of is also invariant under $ as 
follows straightforwardly from the reversibililty of $ : 

/iph^-HM x [0,i/)}) = ft ph ([0,x) x [0,y)) , 

is equivalent to 

£ph([0,x) x [0,y)) = /W${[0,aO x [0,y)}) . (39) 

That is, we may think of /x p h as a repelling measure for $, in the sense that, while it is 
indeed invariant, any slight deviations from this measure, if they are absolutely continuous 
with respect to the Lebesgue measure, will evolve toward the measure /Li p h for the attractor 
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for $ (cf. Eq. (|3lD). In short, we find, for a non-conservative reversible map different 
SRB measures /i p h and /2 p h for the forward and backward time evolutions, respectively. 
And for each time evolution, one plays a role of an attracting measure and the other a 
role of a repelling measure in the sense just explained. This observation is a key element 
of the compatibility between dynamical reversibility and irreversible behavior of statistical 
ensembles. Indeed, when the dynamics is reversible and statistical ensembles irreversibly 
approach a stationary ensemble /i p h for t — > +00, there should exist another stationary 
ensemble /2 p h which is obtained from /i p h by the time reversal operation. However, the 
new stationary ensemble /2 p h is repelling and, thus, its existence is compatible with the 
irreversible behavior of statistical ensembles. 

Further, the measure /i p h is mixing with respect to the forward time evolution $' (t = 
0, 1, • • •). Then it is interesting to investigate the relation between Lyapunov exponents and 
the Kolmogorov- Sinai (KS) entropy for ($,/i p h). The Lyapunov exponents for ($,/2 p h) are 
easily found to be 

K{®,fi P h) = -rlnl-llnr , (40) 

and 

\ y (®,p, ph ) = rlnr + l\nl . (41) 

It is useful to note that the positive (negative) Lyapunov exponent of /i p h can be found by 
changing the sign of the negative (positive) Lyapunov exponent of /i p h- This fact is widely 
used in the analysis of the Lyapunov spectrum of thermostated many particle systems 
In our case this result follows from the observations that 

As($,/Xp h ) = -A^/vO , (42) 

and 

A y ($,/2 ph ) = -A x ($,/x ph ) . (43) 



21 



Thus the Lyapunov spectrum changes sign under the exchange of /i p h and /x p h. Now we turn 
to a calculation of the KS-entropy. First we note that the KS-entropy of ($, yU p h) is 

h K s(®,Vph) = \ x {$,fi P h) , (44) 



which follows from the Pesin's identity |IS |, since /i p h is the SRB measure for the map $, and 
which can also be computed directly from the entropy of the generating partition formed by 
the two elements < x < I and I < x < 1. From the same partition, it is readily seen that 
the entropy of ($, /2 p h) is also 

h KS (<$>, /I P h) = /Vi) • (45) 

Then the difference between the positive Lyapunov exponent and the KS-entropy for ($, /2 ph ) 
is 

K(®, ftph) ~ h KS (®, /2 P h) = -A 2/ ($, /i P h) - A x ($, ^ph) > , (46) 

which is strictly positive for I ^ 1/2 since the right-hand side is just the phase space 
contraction rate. Therefore, the mixing system ($, /i p h) violates Pesin's identity, but satisfies 
Ruelle's inequality ]T5| , as it should be since the measure /2 p h is not an SRB measure for <E>. 



IV. SRB MEASURE FOR A BAKER MAP WITH A CANTOR-LIKE 

INVARIANT SET 

Transport properties are also studied for open Hamiltonian systems with a flux bound- 
ary condition [|I3|-|T5| or with an absorbing boundary condition ||i2| , p0| -[2"6|j . Breymann, Tel 



and Vollmer used an open non-conservative system to study an interrelation between the 
thermostated systems approach and the open systems approach ] . Further, a model used 



by Kaufmann, Lustfeld, Nemeth and Szepfalusky |59[ to investigate deterministic transient 
diffusion is a non-conservative open system. One of the important features of those open 
systems is the existence of escape. So we investigate how escape affects the physical measure, 
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(47) 



by using a Baker map with a Cantor-like invariant set |10,3~C|5"8||. The map is defined on the 
unit square [0, l] 2 (Fig. [?[): 

(x/Z,Aiy), xg[0,Z] 
((x — a)/r, A 2 y + b) , xE[a,a + r] 
where < I < a, < r < 1 — a, < Ai < b and < A 2 < 1 — b. Here we introduce an 
escape for points x G (/, a) U (a + r, 1] and an inhomogeneity. For Ai < I and A 2 > r, or for 
Ai > Z and A 2 < r, the map ^ is partially attractive and partially repelling as discussed in 
the previous section and, for Ai = r, A 2 = Z, and I + r — 1, it is reduced to the previous 
model. It is everywhere attractive for Ai < / and A 2 < r, is conservative for Ai = I and 
A 2 = r and is everywhere repelling for Ai > I and A 2 > r. Note that the last case is possible 
only when there exists an escape: I + r < 1. 

We show that when the initial distribution function p is continuously differentiable in 
x, and a dynamical function (p is continuously differentiable in y and continuous in x, an 
expectation value of the dynamical function <p(x,y) with respect to p t decays exponentially: 



(p(x,y)p t (x,y)dxdy = is / <p(x,y) dxdG(y) p (x, y) dH(x)dy + 0(tA ) , (48) 

[0,1] 2 -^[0,l] 2 "'[0,l] 2 

where the decay rate is equal to the remainder volume per iteration: v = I + r, X' = 
max(/,r)(< v). The functions G and H are (possibly) singular functions defined as the 
unique solutions of de Rham equations : 

-M- 



and 



G(y) 



I 



I + r 



H(x) 



1 . 



I TT X 

ttt h {t 
I 

T+r ' 
r 

r+7 

1 . 



b )+ I 



H 



x — a 



+ 



/ + r 



ye [o,AJ 
y e (Ai,6) 
y G [6, A 2 + b] 
y e (A 2 + 6, l] 

x G [0, Z] 
x G (Z, a) 
x G [a, a + r] 
x G (a + r, 1] 



(49) 



(50) 
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The derivation of Eq. (f48|) is outlined in Appendix B. Now we investigate the implications 
of Eq. P8| ) in cases without and with escape separately. 



A. Baker map without escape 

Consider first the case where there is no escape (I + r = 1). The unit square is then 
contracted onto a set which is a direct product of a Cantor set in the y direction and the unit 
interval, [0, 1], in the x direction. This direct product set is nothing but the strange attractor 
of whose Hausdorff dimension 1 < Hp < 2. In this case, H(x) = x and fl48|) reduces 



to an expression similar to (|3q) for the previous example. And the function G defines an 



invariant mixing measure p p h] 

/i ph ([0,x) x [0,y))=xG(y) . (51) 

The measure /i p h is the one which provides long-term averages of dynamical functions and 
thus, is a physical measure. Indeed, when / p (x,y)dxdy = 1, (H8|) gives 



(v) t = <f(x,y)p t (x,y)dxdy = (f(x,y) dxdG{y) + (t{A'}*) • (52) 

J[0,1] 2 "'[0,l] 2 v ' 

Since /i ph is smooth (more precisely absolutely continuous with respect to the Lebesgue 
measure) along the expanding coordinate x, it is an SRB measure. As shown in Fig. || the 
graph of G(y) is a typical devil's staircase and the measure /x p h is singular. Note that the 
support of ^i p h is the strange attractor of 

B. Baker map with escape 

When I + r < 1, almost all the points escape the unit square and there appears a fractal 
repeller, which is singular both in expanding and contracting directions. Then, the invariant 
measure supported by the fractal repeller is different from the physical measure defined by 
Eq. (|3p, which is not invariant but conditionally invariant under ^ [ f23| . [2""l| 1 . 
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First we consider the physical measure. The formula (f48|) holds even when l + r < 1. By 
setting (p = 1 in (fJHJ), the renormalization factor / p t (x, y)dxdy is found to be 



Pt(x, y)dxdy = v l \ Po(x, y) dH(x)dy + OltX' 1 ) 

[0,1] 2 "'[O,!] 2 



Hence, the expectation value of p at time t is 



(f(x,y)p t (x,y)dxdy / ✓ 

(^)t = ^7 = / n ^ll)«M?(y) + oM- . (53) 

' Pt (x,y)dxdy W V I "J 1 



'[0,1] 2 

This implies that the physical measure defined by Eq. (gj) is identical to that for the Baker 
map without escape: 

p ph ([0,x) x [0,y))=xG(y) , (54) 

provided that the ratio l/r is common in two cases. As in the previous case, the physical 
measure yU p h is singular and is supported by a direct product of a Cantor set C\ along y and 
the unit interval along x: [0, 1] x C\, which is the unstable manifold of the fractal repeller. 
Since it is smooth along the expanding x-direction, it is an SRB-like measure (in the sense 
that, although not invariant, it is smooth along the expanding direction). 

We observe that the support of /i p h is the unstable manifold of the repeller and is not an 
invariant set with respect to Accordingly, the physical measure p p h is not an invariant 
measure. Indeed, the measure p p h satisfies 

^(rHlO.x) x [0,i/)}) = (l + r) /i ph (j0,x) x [0,y)) , (55) 

which implies that /i p h shrinks as time goes on. This can be seen immediately from the 
functional equation for G. For example, when y e [0, Ax], 

/iphfVHM x [0,y)}) = fi ph ([0,lx) x [0,y/Ai)) = lxG{y/^) 

= (l + r)xG{y) = {l + r) p ph {[0, x) x [0, y)) . 

Then, since the measure p p ^ satisfies 
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for any Borel set E C [0, l) 2 , it is conditionally invariant p3| , p4 



Now we turn to an invariant measure yUj n on the repeller, which is defined by 

/i in ([0,x) x [0,y)) = H(x)G(y) . (56) 

The invariance can be seen straightforwardly from the functional equations for G and H. 
For example, when y G [b, A 2 + b] , one has 

/^(^{[O, x) x [0, y)}) = /i in ([0, Ix) x [0, 1) U [a, rx + a) x [0, (y - 6)/A 2 )) 

= /i in ([0, Zsc) x [0, 1)) + // in ([a, rx + a) x [0, (y - 6)/A 2 )) 

= H(lx)G(l) + {H(rx + a) - H(a)}G((y - 6)/A 2 )) 

= H(i)r!- + J ^ r H{x)G{{y - 6)/A 2 )) 
/ + r Z + r 

= #(£)G(y) = /i in ([0,x)x[0,y)). 

Since is a fractal function similar to those discussed in Sec. II, the invariant measure is 
singular both along the contracting and expanding directions. This can be easily understood 
as follows: since the map ^ eventually transforms the unit square into the unstable manifold 
of the repeller, which is a direct product of a Cantor set C\ along y and the unit interval 
along x: [0, 1] x C\, the measure becomes singular along y. On the other hand, only the 
orbits starting from the stable manifold of the repeller remain in the unit square and the 
stable manifold is a direct product set of the unit interval along y and a Cantor set C 2 along 
x: C 2 x [0, 1]. Thus, the invariant set is a subset of C 2 x [0, 1]. As a result, the invariant 
measure becomes singular also along x. Actually, the direct product of the two Cantor sets 
C 2 x C\ is the fractal repeller of \& '. 

In a similar argument to the derivation of (f48|) , one can show that the invariant measure 
/i; n is mixing with respect to ^l. As shown in Appendix C, the invariant measure [i m is a 
Gibbs measure HE!1,|2H2| • 
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V. A BAKER-TYPE MAP UNDER A FLUX BOUNDARY CONDITION 



To illustrate a stationary state for an open system with a flux boundary condition, we 
study a simple model of a chemical reaction. In simple reactions such as R <-> I P, the 
reactant R, the intermediate I and the product P consist of the same atoms, and they can 
be specified by configurations of atoms, or points in the atomic-configuration space. An 
example of the reaction R <-» I <-> P is an isomerization, or a change in the conformation 
of a molecule, such as the "chair" to "boat" transformation of cyclohexane, where R is 
the chair-shaped isomer, P is the boat-shaped isomer and I is an intermediate unstable 
isomer, all of which consist of six carbon atoms and twelve hydrogen atoms. Hence, the 
reaction process can be regarded as a dynamical process where each trajectory starting 
from a reactant state to a product state represents an individual reaction R —>■ P (cf. Ref. 



60| and references therein). A connection between chemical reactions and the escape-rate 



formalism was investigated by Dorfman and Gaspard p2| and a Baker-type model of a 
chemical reaction was studied by Elskens and Kapral |40| . In this section, we introduce a 
Baker-type model of a reaction R «-> / «-> P and study its statistical properties for two 
cases: In one case, the system is closed, and in the other case, a flux boundary condition is 
imposed. Now we begin with the phenomenology. 



A. Phenomenology 

For a chemical reaction R <->• / <-> P, the discrete-time version of the phenomenological 
rate equation is 

Rt+i = k RR R t + km I t , 

It+i = k IR R t + k n I t + k IP P t , (57) 
Pt+i = kpil t + kppP t , 

where Rt, It, and Pt are concentrations of the reactant R, intermediate J, and product P, 
respectively, and kAB (A, B = R,I, or P) are rate coefficients. Since the sum R t + It + Pt 
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is preserved, the rate coefficients satisfy a sum rule : 



kpp + k IR — 1 , 
km + k H + k PI = 1 , (5£ 
kjp -\- kpp = 1 . 



The stationary state solution of Eq. (p7|) is then given by 

-Rst fciy P s t fcp/ 



/ / ' / /• ■ (59) 

1 st -<st MP 



Now we consider a stationary solution of Eq. (p7|) under a flux boundary condition, 
where the concentrations of the reactant R and the product P are fixed to given values 
R ex and P ex , respectively. This may be realized e.g., by introducing source terms to the 
equations for the reactant and product and by adjusting them so as to keep the values of 
R t and P t constant. Eq. (|57D is then reduced to 



It+i = k IR R ex + k H I t + A;/pPex , (60) 
which has a stationary solution 

k IR R ex + kipP cx . . 

^fl = ; ; • 61 

1 - k n 

The deviation 5I t = It — la of the intermediate concentration from the stationary value 
decays exponentially : 

SI t = k n Sit-! — • ■ • — kjj SI ■ (62) 

The stationary state admits a non-vanishing concentration flow from the reactant to the 
product : 

t it t,r> kpikiR k RI k IP 

JR-^P — Kpilfl — ^IP-Lex — Z } -Kcx _ " ; -Tex , 

1 - k H 1 - k n 

and, hence, can be regarded as a stationary state under a flux boundary condition. 
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B. A Baker-type model and a stationary state for a closed system 



We introduce a Baker-type model of the reaction process R <-> I <-> P. Microscopic 
dynamical states of each species R, I, or P are represented by the points in a unit square 
(0, l] 2 . Area preserving asymmetric Baker maps are used to describe the dynamics of the 
reactant and product states, and a Baker map similar to the one discussed in the previous 
section is used for the intermediate state dynamics. The model is then defined by (Fig. |]): 
' (l:f,A lV ) , xe(0,l) 
^:|5| > ( 6 _A 1 )^ , xe(l,a] 

(i : ^=-^,A 2 y + fc) , xe{a,a + r] 



*'(/ : x,y) 



f Zgiy ,(l-6-A a )y) , xe(o + r,l] 



(64) 



(^:yzr^,(6-Ai)y + Ai) , are(0,6-A a ] 
(^=f5|^ > (i-& + Ai)y + 6-A 1 ) , xe(b-A 1 ,l] 



(65) 



*'(P:^2/) 



J : i _ b- A > (1 - & - A 2 )y + 6 + A 2 ) , x e (0, 1 - 6 - A 2 



(66) 



( P: X ~ 6 1 |^ 2 +A2 ,(& + A 2 )y + l-6-A 2 ) , ze(l-&-A 2 ,l] 

where (x, y) denotes a point in a unit square and a (a — R, I, or P) distinguishes different 
species. 

As before, we study the time evolution of a measure starting from an initial measure 
which is absolutely continuous with respect to the two-dimensional Lebesgue measure. The 
Frobenius-Perron equation for the density function p t (a : x, y) (a = R,I, or P) is obtained 
from 



£ / dxdyS ((a : x , y) - *'(/? : y')) M : x> y>) , 
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where the delta function 8 ((a : x, y) — (ft : x', y')) stands for the product 5 a ,p8(x — x')8(y — 
y'). By integrating the Frobenius-Perron equation with respect to y, one obtains the evo- 
lution equation for the partially integrated distribution function Qt(ot : x,y) — dy'p t (a : 



x,y') 



Qt+i(a :x,y) = FQ t (at : x, y) + R t (a ■ x, y) , (67) 

where a contraction mapping T and a functional Rt of Qt(oi '■ x, 1) are given in Appendix 
D (cf. Eqs. 0, (p3|), (PI ) and Eqs. flPgj) , (|D6|), (|D7|), respectively). 



Now we consider the case where the initial measure is uniform along the expanding x 
direction and thus, Qo(a : x, y) does not depend on x. Then, as seen from the expressions 
of T and R\, the partially integrated function Q t (ot '■ x,y) at time t is also independent of 
x. Particularly, its value Q t (a : x,y — 1) = Qt(ot) at y — 1 obeys 

Q t+ i(R) = (1 - b + Ax)Q t (R) + (a- l)Q t (I) , 

Q t+x (I) = (b- Ai)Q t (#) + (/ + r)Q t (/) + (1 - b - A 2 )Q t (P) , (68) 
Q t+1 (P) = (1 - a - r)Q t (I) + (b + A 2 )Q t (P) . 

Since the distribution is uniform along the x direction, Qt(ce) is equal to the total probability 
of finding the system in a species a : Qt(az) = J dxdyp t (a : x, y). Therefore, Eq. (|68f) has 
exactly the same form as Eq. (|57D where the corresponding rate coefficients are given by 

kRR = 1 - b + A x , k RI = a-l , 

km = b - A x , k H = l + r , k IP = 1 - 6 - A 2 , (69) 

fcp/ = 1 — a — r , fcpp = 6 + A 2 . 



Note that the rate coefficients given above satisfy the sum rule Eq. ( |58"D and hence, the sum 
Qt{R) + Qt{P) + Qt(P) is constant in time. We also remark that the rate coefficients ( |6"9~D 
admit a simple geometrical interpretation. As an example, we consider which is the 
transition probability from the intermediate I to the reactant R. According to the definition 
of the map a rectangle (I, a] x (0, 1] moves from the intermediate states to the reactant 
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states and, thus, its Lebesgue area (a — I) corresponds to the transition probability kju from 
I to R, or we have the second equation of (|69|). The other rate coefficients can be obtained 
in the same way. 

According to the discussion given in the previous subsection, the stationary solution of 
1) ^ 



Q*t(R) 



a — I 



Q s t(l) , Q s t(P) 



1 — a 



b- A 1 

and the corresponding distribution is given by 

a — I 



1 - b - A. 



-Qst(/) 



(70) 



Q st (R : x,y) 



-Q st (I)q st (R:y) 



b- Ai 

Q st (P :x,y)= IZ^ll Qst(I)q s t(P ■ y) 



Q st (I :x,y) = Q B t(I)q s t(I ■ y) , 



(71) 



where the singular functions g s t(« : y) (« = R,I, and P) are the unique solutions of de 
Rham equations 



Q s t(I ■ y) 



I U {1:1) 
l + (a-l)q st (R:^E^ 



a + r q st [I : 



"Al 



a + r+(l-a- r)q st (p : f _ 5 _ ^ 



ye (o,Ai] 
ye (Ai,6] 
y e (6, A 2 + b] 
ye(A 2 + b, 1] 



(72) 



(&-A0g st J: 



6- Ai 

6 _ Al + (1 _ & + AOfct (i2 : f Ift + Ai 



ye (0,6- Ax] 
j/e(6-A!,l] 



(73) 



?st(P : 2/) 



■I- - ^ - A 2 + (h -r A L .)(/ S| (P; ^l|±|±A2 



(74) 



;i - 6 - A 2 )g st [1 : 1 _^_ A J , y e (0, 1 - b - A 2 ] 

ye(l-6-A 3) l] 

One can show, by exactly the same method as before, that the partially integrated function 
Qt(ot : x, y) approaches the stationary state Q s t{a : x,y) provided that the initial function 
Qo(a : x, y) is continuously different iable with respect to x. Then since the asymptotic 



31 



measure Q s t(at : x, y) is absolutely continuous with respect to the Lebesgue measure along 
the expanding x direction, it it the SRB measure. Note that, in this case, the quantity 
Q s t(/) is a functional of the initial distribution Qo(a : x, y). 



C. A stationary state under a flux boundary condition 

As shown in Sec. VA, one has a flux from the reactant to the product when their 
concentrations are fixed to the values different from the equilibrium ones. This can be 
realized in the Baker-type model by fixing the measures of the unit squares corresponding 
to the reactant and the product to uniform Lebesgue measures with different densities. So 
we set Qt(R : x,y) — R cx y and Qt(P '■ %, y) — P cx y (t = 0, 1, •■•). Note that the same 
procedure was used to achieve the flux boundary condition for the finite multi-Baker chain 
||13|| . Then, we have only one time-dependent variable Qt{I '■ x,y), which will be abbreviated 
as Q t (I : x, y) = Q t (x,y). Then the equation of motion ( p7| ) for the partially integrated 
distribution function Q t becomes 

Q t+1 {x,y) = rQ t {x,y) + R t {x,y) + S{y) , (75) 

where the contraction mapping T' and Qt(x, Independent part R t (x,y) are given by Eqs. 
(|B3|) and (|B4|) of Appendix B, respectively, and the source term S(y) is due to the reactant 
and product states 

fo, ye[o,Ai] 
R ex (y-A 1 ), ye(A u b) 



S(y) 



(76) 

R^ib-Ax), ye[b,A 2 + b] 



, R ex (b - Ai) + P cx (y - b - A 2 ) . y E (A 2 + b, 1] 
By exactly the same argument as before, Eqs. ([75]) and (|76| ) is found to have a solution 



Qt(x,y) = b Al R cx r) R (y) + 1 b A2 P cx Vp(y) 
1 — I — r 1 — / — r 



W f Q dH(x'){Q (x', 1) - 1 & _ / Al r R cx - \ ^ _ A " P cx ]G(y) + 0(tA rt ) , (77) 
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where G(y) and H{x) are singular functions introduced in Sec. IV (cf. Eqs. (f49|) and (|50| 
and the functions r] R (y) and rjp(y) are the unique solutions of de Rham equations 

ivn(jL), ye[o,Ai] 

1 - r+ 1 fe~-A 1 r (y~ 6 ) ' ^ e ( A i' 6 ) 

r rjni^^-) + 1 - r , y G [6, A 2 + b] 

l , !/e(A 2 + M] 



(78 



l+ , 1 ~ f -r (y 



r 



ye [o,AJ 
2/ e (Ai,6) 

y G [b, A 2 + b] 
y G (A 2 + 6, 1] 



1-T^Y 2 

Therefore the stationary measure /ifj for this open system is given by 

b- A 2 



/i fl [0,x) X [0,y) 



b ~ A l n / \ 1 

■ R ex x r] R (y) + — 

L — / — r 1 



Ax a; r] P (y) 



(79) 



(80) 



As it is smooth (i.e., absolutely continuous with respect to the Lebesgue measure) along 
the expanding coordinate x, it is an SRB measure. However, the measure /ifi is different 
from the SRB measures obtained before since it is absolutely continuous with respect to 
y. It is only weakly singular in the sense that the density is discontinuous on a Cantor- 
like set (cf. Fig. ITUp. Actually, such stationary measures become truely singular only for 
infinite systems It is interesting to observe that the conditionally invariant measure 

A*phf[0) x ) x [0) y) \ — x G(y) appears in the time evolution of the measure (|77|) as the decay 
mode, i.e., as the Pollicott-Ruelle resonance. 

Now we investigate macroscopic aspects. The total measure /i t nO, l] 2 ^ = Jq 1 dxQ t (x, 1) 
correponding to the concentration of the intermediate I is 

to f(0, l] 2 ) = *~ Al _ Re, + \~ h 7 k ! Ax + SSflo + 0(tA'*) 



/ — r 



/ — r 



where 5/j.q stands for the deviation of the initial measure from the stationary state /xg 
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5fi = C dH(x')\Qo(x', 1) - * , Al R cx - \ \ Aa P e: 
Jo I 1 — I — r 1 — I — r 

The stationary state /Xfj admits a flux Jr^p from the reactant state to the product state. 
Indeed, the measure /ifl((a + r, 1] x (0, l]j is transfered to the product state and the measure 
(1 — b — A2)-P ex is transfered from there. This implies the existence of the flux Jr^p given 
by 



J r ^p = /i fl {{a + r, 1] x (0, 1]J - (1 - b - A 2 )P CX 

(l-a-rXfe-Aj) ( a -f)(i_ 6 _A 2 ) 

1 - Z - r ex 1 - / - r ex ' [ } 

With the specification fl39f) of the rate coefficients, the expressions of the stationary 



measure and the decay mode as well as (|82D of the flux agree with their phenomenolog- 



ical counterparts (|6T|) , (|62]) and (|63]), respectively. In short, the measure /ifi describes a 
nonequlibrium stationary state which has a non- vanishing flux Jr^r. 



VI. CONCLUSIONS 

We have explicitly constructed SRB measures for three Baker-type maps employing the 
similarlity between the de Rham equation and the evolution equation of partially integrated 
distribution functions. Now we give a few more remarks. 

a) In our examples discussed in Sees. Ill and IV, we have encountered different types 
of contracting dynamics that lead to invariant states all of which have singular measures 
supported on fractal sets. In the first case we considered a non-conservative reversible 
Baker map on the unit square which globally preserves the area of the square but forms 
an attractor-repeller pair due to the local contraction and expansion properties of the map. 
This attractor is fractal since it has the information dimension less than 2, but it is dense in 
the unit square. In the second example of the Baker map with a Cantor-like invariant set, 
we considered two different cases - global contraction onto a fractal set with and without 
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escape of points from the unit square. When there is no escape, the invariant set is a fractal 
attractor which is a non-dense subset of the unit square. If we add the possibility of escape, 
then the invariant set is fractal in both x and y directions. 



b) It is worth mentioning that one can convert a physical measure for a system with escape 
to the proper invariant measure on the repeller by incorporating into the averaging process 
described in Eq. (|]) both the characteristic function on the region of the phase space from 
which points escape, as well as a "survival" function which is unity if the phase point is in 
the region for the time interval < r < T, with T > t, and zero otherwise. Then by taking 
the limit T —>■ oo, one recovers the invariant measure on the repeller. Such a procedure was 
used by van Beijeren and Dorfman in order to correctly compute the Lyapunov exponents 
on the repeller for a Lorentz gas [p^j . This procedure is closely related to one described by 



Hunt, Ott, and Yorke [EH] to obtain natural measures on invariant sets. 



c) We have encountered three different physical measures, which are characterized by the 
smoothness along the expanding direction. For closed non-conservative systems, the physical 
measures are singular and invariant. For an open system under an absorbing boundary 
condition, the physical measure is singular and conditionally invariant. And for an open 
system under a flux boundary condition, the physical measure is invariant, but absolutely 
continuous with respect to the two-dimensional Lebesgue measure. It is only weakly singular 
in the sense that the density is discontinuous on a Cantor-like set. 

It is interesting to note that the conditionally invariant measures appear as the decay 
modes (i.e., the Pollicott-Ruelle resonances) for an open system under a flux boundary 
condition. Such a relation exists more generally. The Pollicott-Ruelle resonances are defined 
as the functionals T acting on a dynamical variable <p |n,3"2|-3"£|, which satisfy, in case of a 
map, 

(r,^o$*) = c* (i» , 
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where $ is a map and ( < 1) is a decay rate. When the characteristic function xa of 
any Borel set A is in the domain of the functional T and (r, xm) with M the whole 
phase space, then the (possibly complex) measure defined by 



n{A) = (r, XA ), 

is conditionally invariant since yu($ _1 yl)//i($ _1 M) = fi(A). Such examples are the hydro- 
dynamic modes for open Hamiltonian systems p^ , |36|j37[| . 

d) For a non-conservative reversible map we find different SRB measures /i p h and /Z p h 
for the forward and backward time evolutions, respectively. For each time evolution, one 
plays a role of an attracting measure and the other of a repelling measure. This is a key 
element of the compatibility between dynamical reversibility and irreversible behavior of 
statistical ensembles. Indeed, when the dynamics is reversible and statistical ensembles 
irreversibly approach a stationary ensemble /i p h for t — > +00, there should exist another 
stationary ensemble /x p h which is obtained from /i p h by the time reversal operation. Since 
the new stationary ensemble /x p h is repelling, its existence is compatible with the irreversible 
behavior of statistical ensembles. A similar situation was observed for open Hamiltonian 



systems under a flux boundary condition [13,14|. Such systems admit two different stationary 



states, one is the time reversed state of the other. In this case, an attracting stationary state 
for the forward time evolution obeys Fick's law and a repelling state obeys anti-Fick's law. 



e) As mentioned in Introduction, an SRB measure for a map may be constructed as follows 
pOj : 1) Approximate the measure by iterating an initial constant density finite times, 2) 
calculate the average with its result and 3) take the limit of infinite iterations. On the other 
hand, in our construction of an SRB measure, a functional equation for the partially inte- 
grated distribution function similar to the de Rham equation is derived from the evolution 
equation for measures. Note that the functional equation is a direct consequence of the self- 
similarity of the measure. An iterative method of solving the functional equation is similar 
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to the procedure explained above. However, our method has some advantages. Firstly, the 
de Rham-type functional equations are fixed point equations for contraction mappings and, 
as a result, the iterative solution strongly converges to the limit. On the contrary, the iter- 
ative approximation of the density function does not converge by itself. Therefore, one can 
obtain a good approximation to the cumulative distribution function of the SRB measure 
by less iterations than by the conventional method. Secondly, in the functional equation 
approach, one can systematically extract exponentially decaying terms which are typically 
higher order derivatives of singular functions in the sense of Schwartz' distributions. Thirdly, 
although the measure is defined as the solution of a functional equation, the average values 
of certain dynamical functions can be calculated analytically as illustrated in Sec. II. So 
far, the functional equation method was mainly applied to piecewise linear one- dimensional 
maps and Baker-type maps |13| , |37| , |38| , but we believe that the method can also be applied 
to other systems if the expanding and contracting directions can be well-separated. 
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APPENDIX A: PROOFS OF EQ. Q AND THE LEMMA 



1. Derivation of Eq. (|19|) 



As explained in Sec. Ill, the evolution equation fllgp for the partially averaged distribu- 
tion function Pt(x, 1) = Jq dy'p t (x,y') 

P t+1 {x, 1) = r P t {r x + l,l) + l P t (l x, 1) , (Al) 

is the Frobenius- Perron equation for a one- dimensional exact map Sx = x/l (for x G [0,/]) 
and Sx = {x — l)/r (for x G which admits the Lebesgue measure as the invariant 

measure. As a consequence, the integral of Pt(x, 1) over the unit interval [0,1] is invariant. 
Indeed, by integrating Eq. (|A1|), we have 



/ dxP t+ i(x, 1) = / dxP t (x, 1) = • • • = / dxP (x, 1) . 
Jo Jo JO 

On the other hand, since po and, hence, Pq is continuously differentiate in x, Eq. ( |A1|) im- 
plies that Pt(x, 1) (t — 1, 2, • • ■) are also continuously differentiate in x. And the derivative 
d x P t (x, 1) satisfies the equation 

d x P t+1 (x, 1) = r 2 d x P t {r x + 1, 1) + I 2 d x P t {l x, 1) , 
which leads to an inequality 



\\d x P t+1 (; 1)|| < (r 2 + I 2 ) \\d x P t {; 1)|| < X\\d x P t (; 1)11 < • • • < X t+1 \\d x P (; 1)|| , 

where the function norm is defined by \\d x P t (-, 1)|| = sup xe [ 01 ] \d x P t (x, 1)| and A 
max(l, r) < 1. Since the function Pt(x, 1) can be represented as 
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P t {x, 1) - C dx'P t (x\ 1) = f 1 dx' x 1 d x ,P t {x' : 1) - f 1 dx'd x ,P t (x', 1) 

JO JO Jx 

we finally obtain 



\Pt(x,l) 



J dx'P (x', 1)| < dx' x' + j dx'}\\d x P t (-, 1) 



<l\\d*Pt(; 1)|| <l\*\\W, 1)||, 



or Eq. dH). 



2. Proof of the lemma 

Since the statement of the lemma given in the text is not technically complete, we first 
give the precise statement and then prove it. 

Lemma Let J 7 be a linear contraction mapping with the contraction constant < A < 1 
defined on a (Banach) space of bounded functions equipped with the supremum norm || • |j 
(i.e., < A||/||j. And let g t = + v t g^ + g[ be a given function where , g^ and 

g t are bounded functions, v is a constant satisfying X < \u\ < 1, and \\gl \\ < KX with 
some constant K > 0. Then, the solution of the functional equation 

ft+i = Fft + 9t, (A2) 

is given by 

ft = £ 0) + Jf£ + 0(f A') , (A3) 
where f{® and f£) are the unique solutions of the following fixed point equations 

/i 0) = ^ + 9 {0) , (A4) 
/£> = i-F/i? + ^ • (A5) 
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The proof of the lemma is straightforward: From Eq. one obtains 



t-i 



s=Q 

By rewriting the right hand side in terms of g^\ g^> and g^\ this relation leads to 

oo oo 

s=0 s=0 

, t— 1 OO OO n 

+{^/o + E ^2L. - E ^V 0) - ^ E v-^rM . 

s=0 s=t s=t J 

The first sum Y^Lo F s g^ in the right-hand side is f$ and the second sum Y^Lo v~ s ~ l F s g^ 
is , which satisfy Eqs. (|A4j) and fl^5| ) respectively. By repeatedly using the property of 
T, we have 



II/* - /i 0) - ^11 < A ll/oll + f^j + 



if t 



which is 0(iA*) and implies the desired result 



APPENDIX B: PHYSICAL MEASURE FOR A BAKER MAP WITH 
A CANTOR-LIKE INVARIANT SET 



In this Appendix, we outline the derivation of Eq. fl48|), which is quite similar to that of 



Eq.( 



^From the definition (47) of the map one finds that the Frobenius-Perron equation 
governing the time evolution of the distribution function p t (x,y) is given by 

' aT^'It) ' ye[o,A!] 
o, ye{A l7 b) 

-£-p t (rx + a,y-^-) , ye[b, A 2 + b] 
10, y E (A 2 + 6, 1] 

which leads to a contractive time evolution of the partially integrated distribution function 
Qt(x,y) = $ dy'p t (x,y'): 



p t+ i(x,y) = < 



(Bl) 
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Q t+1 (x,y) = F'Q t {x,y) + Rt{x,y) 



(B2) 



where T 1 stands for a contraction mapping: 



lQt(lx, j . 



F'Qt(x,y) 



rQ t {rx + a, 

Lo, 



and Rt(x,y) is a function of Qt(x, 1): 

ro, 



Rt(x,y) 



lQt(lx,l) 
lQt{lx,l) 



ye [0,A X ] 
2/ e (Ai,6) 
y e [6, A 2 + 6] 
2/G(A 2 + 6, 1] 

ye [o,AJ 
y e (A Xj 6) 
y G [6, A 2 + 6] 



(B3) 



(B4) 



as 



t ZQi(/x, 1) + rQ t (rx + a, 1) . y G (A 2 + 6, 1] 
In terms of Qt(x,y), the expectation value of a dynamical variable tp(x,y) is expressed 



(B5) 



(p(x,y)p t (x,y)dxdy = / <p(x,y) dxd y Q t (x,y) 

[0,1] 2 J [0,1] 2 



where d y stands for the Riemann-Stieltjes integral of Q t with respect to y |53 



To solve Eqs. ( |B2| ), ( |B~3| ) and (P4|), we first investigate the equation of motion of Qt(x, 1): 



1) = lQ t (lx, 1) + rQ 4 (ra; + a, 1) . 



(B6) 



We observe that, in terms of a function H(x) defined by a de Rham equation (|50D, one has 



dH{x')Q t {x', l) = u f 1 dH{x')Q t ^{x', 1) • • • = z/* T dH{x')Q {x', 1) , (B7) 

JO JO Jo 

Qt{x,l)- I dH{x')Q t {x',l) = £dx'(H{x') - 6{x' - x)\d x ,Q t {x' ,1) , (B8) 



where v = I + r(< 1) and A' = max(7, r)(< 1). On the other hand, Eq. 
inequality 



leads to an 



RQtO, 1)11 < A'H^Qt-iC-, 1)|| < • • < X'^Qoi; 1)|| , 
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with \\d x Q t (', 1)|| = sup^g^y \d x Q t {x, 1 



) | . Then, one obtains from Eqs. ( |B7| ) and (p8| ) 



Q t (x, 1) = i/ [ dH(x')Q (x', 1) + 0(A") 



and, thus, 



Rt{x,y) 



v 



/ / dH(x')Q (x',l)g^(y) + O(\ ft ) , 
Jo 

fO, ye[0,A 1 ] 



(B9) 



{1) (y) 



< l , y G (Ai,A 2 + 6] 
J + r . y G (A 2 + b, 1] 



(BIO) 



Since v = I + r > max(/, r) = A', the lemma of Sec. Ill can be applied to the evolution 
equations ([B~2D, flB"5| ) and ( |B~4] ) and one obtains 



where G(y) is the solution of a de Rham equation fl4~9|) . The desired result (0) immediately 
follows from Eq. (P 1 1|) . 

APPENDIX C: ON AN INVARIANT MEASURE /x in FOR A BAKER MAP 

WITH ESCAPE 

In this appendix, we show that the invariant measure /i; n on the fractal repeller considered 
in Sec. IV is a Gibbs measure. 

To show that, we first observe that the image \l/ m [0, l] 2 of the unit square by the map 
ty m consists of 2 m horizontal strips, which will be referred to as if| m ' , , • • • H^m ; and 
that the pre-image of \l/ _ri [0, l] 2 of the unit square by the map ty n consists of 2 n vertical 
strips, which will be referred to as V^, V^, ■ ■ ■ V^m. Then, boxes ifj m ' fl Vj n ^ generated by 
the overlap procedure provide a generating partition of the repeller. As easily seen, for each 
box i/] m ' H Vj , a stretching factor for a time interval [—m, n — 1] 





t=—m 



(CI) 
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does not depend on the initial point (x, y) G H^nVj , where A x (x, y) is the local expanding 
rate defined by X x (x, y) = — InZ (if x G [0,/]) and X x (x,y) = — lnr (if x G [a, a + r]). Note 
that, when (x, y) G ffj m -' fl V^"', the pre-image \l/ _fc (x, y) is unique for k = 1,2, • • -m. We 
show that the /ii n -measure of the box if j m ' fl v} n ^ is given by 



p—Uij (n,m) 

fi in (Hl m] n Vf ] ) = ^ , (C2) 

' * J ' p—Uij(n,m) ' v 7 

which implies that the measure is a Gibbs measure |T|,|H|iT),2"4]]. Note that, since the 
numerator of Eq. (|C2|) is a product l s r n + m - s Q f (n + m) factors (s = 0, 1, 2, • • -n + m) and 
there are (n + m)!/{s!(n + m — s)!} boxes with this value, the sum of the numerators, or 
the normalization factor, is 



n+m / , \| 

„— Uij(n,m) _ \ " V ~ r is n+m— s (i _i_ \ra+m _ -(n+m)fc 



where n = — ln(/ + r) is the escape rate 

Before proving Eq. ( p2|) , we verify it for a simple case n = 2 and m = 1. In this 
case, we have V? 1 = [0, Z 2 ] x [0,1], vj 2] = [a,rl + a] x [0,1], vf ] = [/a, Zr + /a] x [0,1], 
Vf ] = [ra + a,r 2 + ra + a] x [0, 1]; = [0, 1] x [0, Ai] and = [0, 1] x [6, A 2 + 6]. 
As an example, we consider a box fl vj 2 ^ = [la,lr + la] x [b, A2 + 6]. For any point 
(x, y) G i?2 n W 2 ') we nave < x < Z and \I/(x, y), 1 i r_1 (x, y) G [a, a + r] x [0, 1] and, thus, 

u 2)3 (2, 1) = J] A, (V(x, y)) = -2 In r - In I . 

t=-i ^ ' 

On the other hand, from the functional equations of H(x) and G(x), we have 

/^(tff 1 H Vf ) = {#(Zr + la) - H(la)}{G(A 2 + b) - G(b)} 

' ' {H(l)-H(0)}-^-{G(l)-G(0)} 



I + r I + r I + r 

lr 2 _ exp(- M2 , 3 (2, 1)) 
" (l + r) 3 (l + r) 3 

which is 
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Now we go to the proof of Eq. (|U2|) . Since one can write = [0, 1] x [aj m ',/3j"^] and 



yW = [ 7 f,$f]x[0,l],and thus, 



^{Ht n] n vf ) = {G((3t ] ) - G{a [ r ] )}{H{5f ] ) - H^f)} 

= ^ m (Hl m] )nUv} n] ) , 



it is enough to show 



and 



r i e ~ u j( n ) 

MVh = , (C3) 



p -Ui(m) 

^(Hl m] ) = jj—^ , (C4) 



where ttj(n) = J2t=o ^ M^O^ y)J and Ui{m) = Ylt=-m ^xi ^(x, y)j are stretching factors 
for a vertical strip and a horizontal strip -ffj m ', respectively. As before, the stretching 
factors are constant on each strip. The relations (|C3|) and ( |C4|) are proved by induction with 
respect to n and m, with the aid of the functional equations for H(x) and G(y), respectively. 
Since the proofs of the two relations are similar, we show ( |C3| ) only. 

It is easy to see that ( |C3| ) holds for n = 1. Now we suppose that Eq. (|C3| ) is valid for n. 
As easily seen, a vertical strip vjj 1 ^ is expressed by some vertical strip Vj as 

v}? +l] = R a n v- Vf 1 

where a = or 1 with R = [0, 1] x [0, 1] and Ri = [a, a + r] x [0, 1]. In terms of 7 j and 
we have 



[1^,16^] x [0,1], a = 

[r 7 j n] + a, r5f ] + a] x [0, 1] . a = 1 
When a = 0, from the functional equation for H(x), one obtains 



T/[n+l] 



/ r , expf— uAn) + In / ) 

rh^ ]) = WrA- • (C5) 
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Then, for (x, y) G Vjj 1 , one has (x, y) G Rq and ^f(x, y) G Vj. Therefore, X x (x, y) — — In/ 
and 

u^n) -In/ = ^A a: (Votf(x,y)) + A x (x, y) = £ A x (V(x, y)) = 1^/(71 + 1) . (C6) 
Similarly, one can verify Eq. (|C6|) when a = 1. Hence, 

n , [n+i] exp(-ry(n + l)) 

/W*' J - (/ + r) n + l 

or Eq. (|C3|) holds for n + 1 and, by induction, it is valid for all positive integer n. 



APPENDIX D: THE EVOLUTION EQUATION OF MEASURES FOR 

A REACTION MODEL 

In this Appendix, we write down the evolution equation of the partially integrated dis- 
tribution function Q t (ot : x,y) — J* dy'p t (a : x,y') (a = R,I, or P) for a chemical reaction 
model introduced in Sec. V. The density function p t (a : x, y) (a = R,I, or P) at time t is 
given by 

p t (a:x,y)= ^ / dxdyS ((a : x, y) - *'*(/3 : x 1 , y')) Po(/3 : x', y') , 

where po is the density function of the initial measure, ^' is the map introduced in Sec. V 
(cf. Eqs. fl64[), (|65f) , and (|66|) ) and the delta function 5 ((a : x, y) — {(3 : x',y')) stands for 
the product 5 a ^5{x — x')5(y — y'). 

By integrating the Frobenius-Perron equation for the density p t (a : x, y), one obtains 
the evolution equation for Q t : 

Qt+i{a : x, y) = TQ t {a : x, y) + #t(a : x, y) , (Dl) 

where a = R,I, or P, a linear contraction mapping T is defined by 
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TQ t {l : x,y) 



lQ t (l:lx,jL) , 

(b - Ai)Q t (R:(b- A l )x, , 



y-b 



rQ t [I : rx + a,—fc 



ye (0,Ai] 
y e (Ai,6] 
y G (6, A 2 + b] 



(l-b-A 2 )Q t (p:(l-b-A 2 )x, y_ b b _^ , ye(A 2 + M] 



(a-Z)Q t (/:(a-Z)a; + Z, r ^ Ar ) , 



TQ t (R:x,y) = { 



ye(0,b-A 1 ] 



(1 - 6 + Ai)Q t (i? : (1 - 6 + Ai)x + 6 - Ai, f _ ^ + ^ ) , J/ e(6-A 1 ,l] 



(l-a-r)Qt ^/:(l-o-r)x + a + r, 1-6 y _ A J , j/£ (0,l-6-A 2 ] 

.PQ t (P::r,y) = < 

(6 + A 2 )Q f (p : (6 + A 2 )x + 1 - b - A 2 , y - I + 1± A2 ) , ye(l-fc-A 2 ,l] 
and R t {oc : x, y) is a functional of <5t(o! : x, 1) : 

fO, ye^AJ 
lQt(I:lx,l) , ye (Ax, 6] 

R t {I :x,y) = { lQ t (/ : Ix, 1) + (b - Ai)Q t (R : (b - A ± )x, 1) , ye (b, A 2 + b] (D5) 
/g t (/ : Ix, 1) + (6 - Ai)Q t (P : (6 - Ai)x, 1) 

+ rQ t (/ : rx + a, 1) , y G (A 2 + b, 1] 



P t (P : x,y) = 



, 

(a - l)Q t (I : (a - Z)x + /, 1) , 



y e (0,6-Ai] 
ye (6-Ax,!] 



(D6) 



, y e (0, 1 - b - A 2 ] 

P t (P:x,y) = (D7) 
1 (1 - a - r)Q t (/ : (1 - a - r)x + a + r, 1) . y G (1 - 6 - A 2 , 1] 

These are the desired results. Note that the contraction constant A of the mapping f is 

given by 
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A = max^/,r, a - 1, 1 - a - r, b - A 1; b + A 2 , 1 - b + Ai, 1 - b - A 2 ^(< 1) . 
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FIGURES 



FIG. 1. Wcierstrass function for a = 2/3 and 6 = 2. 
FIG. 2. Takagi function. 
FIG. 3. Cantor function. 

FIG. 4. Lebesgue's singular function for a = 0.75. 

FIG. 5. Schematic representation of the non-conservative reversible Baker map. The shaded 
rectangle is expanded and the rest is contracted. 

FIG. 6. Partially integrated distribution of the physical measure ^ p h along y for the 
non-conservative reversible map with / = 0.3. 

FIG. 7. Schematic representation of the Baker map with a Cantor-like invariant set. The points 
in the black rectangle are removed. The other shaded rectangles are mapped onto the corresponding 
ones. 

FIG. 8. Partially integrated distribution of the physical measure fi p h along x for the Baker map 
with a Cantor-like invariant set. The parameters are Ai = 0.4, A2 = 0.3, b = 0.5 and l/r = 0.8. 

FIG. 9. Schematic representation of the Baker-type map describing a chemical reaction process 
R <-> / <-> P. Three unit squares describe the dynamical states of the reactant, intermediate and 
product, respectively. The shaded rectangles are mapped onto the corresponding ones. 

FIG. 10. Partially integrated distribution of the stationary measure along x for the reaction 
model under a flux boundary condition. The parameters are chosen as Ai = 0.4, A2 = 0.3, b = 0.5, 
/ = 0.4, r = 0.5, .Rex = 0, and P ex = (1 - I - r)/(l - b - A 2 ). 
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